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ABSTRACT 

The conversion of fast waves to the Alfven mode in a realistic sunspot atmo- 
sphere is studied through three-dimensional numerical simulations. An upward 
propagating fast acoustic wave is excited in the high-/? region of the model. The 
new wave modes generated at the conversion layer are analyzed from the projec- 
tions of the velocity and magnetic field in their characteristic directions, and the 
computation of their wave energy and fluxes. The analysis reveals that the max- 
imum efficiency of the conversion to the slow mode is obtained for inclinations 
of 25 degrees and low azimuths, while the Alfven wave conversions peaks at high 
inclinations and azimuths between 50 and 120 degrees. Downward propagating 
Alfven waves appear at the regions of the sunspot where the orientation of the 
magnetic field is in the direction opposite to the wave propagation, since at these 
locations the Alfven wave couples better with the downgoing fast magnetic wave 
which are reflected due to the gradients of the Alfven speed. The simulations 
shows that the Alfven energy at the chromosphere is comparable to the acoustic 
energy of the slow mode, being even higher at high inclined magnetic fields. 

Subject headings: Sun: oscillations - Sumsunspots - Sun: numerical simulations 

1. Introduction 

The study of solar oscillations has proven to be a powerful tool to infer the properties of 
the solar interior. Global helioseismology, based on the interpretation of the eigenfrequencies 
of the resonant modes of oscillations, has provided a robust description of the internal struc- 
ture and dynamics of the Sun. In the last years, these studies have been complemented by 
focusing on local features. Local helioseismology interprets the full wave field observed at the 
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surface. Several complemetary techn iques have be en developed to probe lo cal perturb ations, 
as Fourier-Hankel spectra l analysis (IBraunlll995r ). ring-diagram analysis (IHilllll988r). time 



dista nce helioseismology ( IDuvall et al. 
1990h - 



19931 ). and acoustic holography (ILindsey fc Braun 



For an accurate interpretation of the measurements obtained from all these procedures, 
it is extremely important to achieve a deep understanding of the physics involved in the 
wave propagation. Due to the increasing interest of helioseismologists in active regions, such 
as sunspots, the knowledge of the inteaction of waves with magnetic structures has been 
greatly developed in the recent years. The conversion from fast-mode high-/3 acoustic waves 
to slow- mode low-/? waves in so l ar active region is well - understood from the theoretical poin t 
of view dCallv fc Bogdanlll993l : Icrouch fc Caiwll2003l : ISchunker fc Callvlbood : ICallvl bood ) . 
An example of the success of this analyt ical development i s the compar ison of the observa- 



tional absorption and phase shift data (B raun et al .11 19881: 



results obtained from the conversion rtCally et al 



200 



3; 



Braun 



Crouch et al 



1995h wi th the modeled 



20051). Numerically. 



the fast-to-slow conversion is also well studied in sunspot-like atmospheres (IBogdan et al 



2003 iRosenthal et al.ll2002l : 



Khomenko fc Collados 2006 



Khomenko et al 



gOOa lFelipe et al. 



201(1) . as well as flux tubes ( Hasan et al. 20031 : Hasan fc Ulmschneider bood ) . See 



Khomenko 



(120091 ) for a review. 



Not all the developments of phenomena associated with wave propagation in magnetic 
structures have reached the same degree of maturity. An insufficiently modelled effect is 
the fast-to-Alfven conversion. Pure Alfven waves can only exist in those mediums which 
are homogeneous in the direction perpendicular to both magnetic field and wavenumber. In 
general, this will not be the situation in a gravitationally stratifed atmosphere with complex 
magnetic field structure, although we will simply refer to Alfven waves even in this case. 
The fast-to-Alfven conversion only occurs when the wave propagation is not contained in the 
same plane of the density stratification and magnetic field and, thus, three dime n siona l (3D) 
analysis are necessary. The study of this process was started by lCrouch fc Callyl (120051 ) . who 
studied the 3D propagation of oscillations in a polytrope permeate d by an uniform m a gnetic 
field of arbitrary inclination, and found downward Alfven waves. ICally fc Goossensl (120081 ) 
obtained that the conversion to the Alfven mode is most efficient for field inclinations from 
vertical between 30 and 40 degrees, and azimuth angles (th e angle between t he m agnetic 
field and wave propagation planes) between 60 and 80 degrees. ICally fc Hansenl (120111 ) found 
that the interaction between fast and Alfven waves is spread across many scale height, unlike 
the fast-to-slow, which is limited around the layer where the sound speed c# and the Alfven 
speed va are similar. As the frequency increases the mode conversion is progressively more 
localized, although at the frequencis relevant to local helioseismology (around 3-5 mHz) the 
fast-to-Alfven conversion region spans the whole chromosphere. 
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Recently, Khomenko and Cally have studied the conversion to the Alfven mode by mean s 
of 2.5D numerical simulations in homogeneo us field configurations (JKhomenko fc Callyll201ll ). 
as well as realistic sunspot-like structures ( Khomenko fc Cally 20121 ). In these works they 
obtained the dependence of the efficiency of the conversion with the inclination and azimuth. 
However, in the sunspot configuration the conversion was only evaluated at some limited an- 
gles correspondi ng to selected 2D p l anes o f the model. The aim of this work is to extend 
the results from iKhomenko fc Callyl (120121 ) to 3D. It will allow us to populate the voids in 
their diagrams of fast-to- Alfven conversion efficiency. More relevant, the development of full 
3D simulations provides the results for the complete wave field, where waves can propagate 
freely in all spatial directions without being restricted to a two-dimensional plane. The con- 
version to the Alfven mode in realistic 3D atmospheres has not been studied before with the 
detailed evaluation of the efficiency discussed in this paper. 



2. Numerical procedures 



The three-dimensional (3D ) non-linear magnetohydrodynamic equations are solved us- 
ing th e numerical code Mancha ( Khomenko fc Collados 2006 : Khomenko et al.ll2008 : Felipe et al 
2010). The code solves the equations for perturbations, obtained after removing the equi- 
librium state from the equatio ns. A Perfect M atched Layer (PML) is used to avoid wave 
reflection at the top boundary (jBerengerlll996l ). while periodic boundary conditions are im- 
posed in the horizontal boundaries. The initial perturbation is set in the bottom boundary 
with small amplitude in order to ensure that the simulations are in the linear regime. 

As a background atmo s phere we use a magnetostatic (MHS) sunspot model, adopted 
from IKhomenko fc Colladosl (120081 ). This model is a thick flux tube with distributed currents, 
azimuthally symmetric and has no twist. We set the height reference z = Mm at the 
photospheric level, where the optical depth at 500 nm is unity in the quiet Sun atmosphere. 
At this height the magnetic field at the axis is 900 G. The spatial resolution is 150 km in 
the horizontal directions and 50 km in the vertical direction. The computational domain 
spans from z = —5 Mm, where the minus sign indicates that it is below the photosphere, 
to z = 2.4 Mm. The upper 500 km (10 grid points) correspond to the PML boundary 
layer, so the effective top of the simulation is located at z — 1.9 Mm. The horizontal 
extent of the domain is x £ [—39, 39] Mm and y £ [—30, 30] Mm, with the axis of the 
sunspot located at x = 0, y = Mm. At 39 Mm from the sunspot axis the thermod ynamic 
variables of the model are taken from Model S (iChristensen-Dalsgaard et al.lll996l ) in the 
deep sub-photosphere layers and VAL-C model (IVernazza et al.lll98lh in the photosph ere 
and chromosphere, stabilized following the method by iParchevsky fc Kosovichevl ( 120071 ) to 
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avoid the convective instability. The axis of the sunspot is given by the lAvrettl (1198ll ) model. 
The atmosphere between the quiet Sun boundary and the umbral model at the axis merges 
smoothly. 

Waves are driven in a few grid points at the bottom boundary at z — —5 Mm. The per- 
turbations in pressure, density, and velocity are calculated analytically as an acoustic-gravity 
wave of a given frequency and way enumber, neglecting the magnetic field and temperature 
gradient (IMihalas fc Mihalasl 11984 ). The detailed form of the perturbations can be found in 
Khomenko fc Callvl h&vK . In this simulation the frequency was set to v = uj/2h = 5 mHz, 
slightly below the maximum cut off frequency reached at the temperature minimum, and 
the horizontal wave number to k T = 1.37 Mm -1 . 



In order to indentify the slow, fast, and Alfven wave modes in the magnetically dom- 
inated region of the computational domain, the velocity and magnetic field perturbations 
have been projected according to their orientation with respect to the equilibrium magnetic 
field onto these three characteristic directions: 
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where 6 is the magnetic field inclination from the vertical and (f) is the field azimuth, mea- 
sured from the x — z plane. The projection e\ ong is along the magnetic field and it selects 
the slow magneto - acoust ic wave in the low-/3 region; the projection e perp was chosen after 



Cally fc Goossensl (120081 ) and it gives the asymptotic polarization direction of the Alfven 
mode; and the last projection e trans is set normal to the other two, corresponding to the fast 
wave in the low-/3 regime. These projections have already been su ccessfully used to separate 



the three wave modes in idealized magnetic fi eld configurations (IKhomenko fc Cally 



as well as more complex magnetic topologies (IFelipe et al.l l2010l : IKhomenko fc Cally 



2011ft 



20121). 



These projections assume that the wavevector k is contained in the xz plane. In these sim- 
ulations k will not be limited to that plane, since the raypath can bend due to background 
variations in the y direction. This would have some effect in the accuracy of the projections, 
but it can probably be neglected. 

As a mea sure of the efficiency of th e conversion to each mode, the time-averaged wave 
energy fluxes (IBray fc Loughheadlll974l ) were calculated in the magnetically dominated re- 
gion. The acoustic energy flux is obtained from the expression: 
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Fac= (PIV), (4) 

while the magnetic energy flux is given by: 

F mag = (Bi x (v x B ))//i . (5) 

where p 1? v, and B x are the Eulerian perturbations in pressure, velocity, and magnetic field, 
respectively, B is the background magnetic field, and /x is the magnetic permeability. In 
the region where va > cs the acoustic energy flux contains the energy of the slow mode, 
while the magnetic flux includes the fast and Alfven modes. The time-averaged energy of 
the three wave modes was also calculated from the relations: 

Elong = P0Cs(vf O ng) ( 6 ) 
E pe r P = PoV A (v 2 perp ) (7) 

Etrans = Po^LrJ ( 8 ) 

where vi ong) v perp) and v trans are the velocity projections into the characteristic directions from 
Equations [TKHl and p is the density in the equilibrium state. These expressions provide an 
approximation of the wave energy assuming equipartition between kinetic and other energies. 
In the case of pure acoustic or Alfven waves, there is strict equipartition between kinetic and 
compressional or magnetic energy, respectively, and Equations ([MED correspond to the real 
energies. 

3. Velocity projections 

When the fast acoustic wave which was driven at the bottom boundary reaches the 
va = cs layer several mode transformations take place. Above that layer, in the magnetically 
dominated atmosphere, the incident wave splits into a slow acoustic mode, a fast magnetic 
mode, and an Alfven mode. 

Figure [I] shows snapshots of the projected velocities scaled with a factor ^poVph, where 
v p h = cs for the V\ ong component and v p h = va for the other two components, after 19 min 
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of simulation. The stationary regime is achieved after about 10 minutes of simulations. The 
inclination of the magnetic field 9 varies from 0° at the center of the sunspot to being almost 
horizontal at the boundaries of the computatial domain. The azimuth = 0° corresponds 
to all the positions at y = Mm with positive X value, and it increases up to = 180° at 
y — Mm and negative X values, including all the angles between these two extremes. 
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Fig. 1. — : Snapshots of the three orthogonal components of the velocity at t — 19 min. The 
blue-red colors mean positive-negative velocity directions; the range of the color coding is the 
same in all panels. The projection v\ ong is scaled with a factor of sj p^cs and the projections 
Vtrans and Vperp with a factor a/poVa- Top panels correspond to horizontal cuts at z = 1.65 
Mm, from left to right: vi ong) v tranS) and v perp . Dashed lines represents contours of equal 
inclination of the background magnetic field, measured at the height where c$ = va- The 
three bottom panels show vertical cuts at y — 7.5 Mm, from top to bottom: vi ong) v transi 
and v perp . Horizontal solid line is the height where cs — va] horizontal dashed line is the fast 
mode reflection level. Magnetic field lines are inclined black lines. 



Each of the waves modes presents a different distribution across the 3D atmosphere 
of the sunspot. In the magnetically dominated atmosphere the slow wave appears in the 
viong projection. The xy cut (top left panel) shows that the conversion to the slow mode is 



-7- 



significative at almost all the positions of the sunspot. However, there are some locations 
where this transformation is specially favoured. The strongest slow wave signal appears for 
X between 3 and 12 Mm and Y between —10 and 10 Mm. This region corresponds to 
moderate inclinations around 9 = 30°. Although the frequency of the wave is below the 
cut-off frequency, the slow mode can reach the upper atmosphere because of the reduced 
cut-off value due to the inclination of the magnetic field. However, the vertical magnetic 
field at the axis of the sunspot avoids the propagation of slow waves around the center of the 
sunspot, and they produce evanescent modes. In general, the amplitude of the slow mode 
in the x-positive region of the sunspot is higher than the x-negative region. The driving 
perturbation generates an acoustic-gravity wave which propagates from left to right. Thus, 
the right half of the sunspot has a better alignment between the direction of propagation 
and the field lines, producing a more efficient conversion from fast acoustic waves (in the 
region where va < cs) to slow acoustic waves (in the region where va > cs)- 
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Fig. 2. — : Variation of the temperature with height at the axis of the sunspot (solid line), 
at 7.5 Mm from the axis (dotted line), and at the quiet Sun atmosphere (dashed line). 

The vertical cut of the slow mode (Figure (TJi) shows an interference pattern above the 
solid line, near the top boundary. It appears mainly in the left part of the sunspot, between 
x = —20 Mm and x = 5 Mm. The interference is produced between the upward propagating 
slow mode and a downward propagating wave produced by the partial reflection of the slow 
mode due to the steep temper ature increase at the chrom osphere. Similar result was found 
in the 2.5D simulations from iKhomenko fc Callyl ( 12012I ). where they confirmed that this 
reflection is a physical effect rather than an artifact from the top boundary. This reflection 
is stronger near the center of the sunspot, where the gradient of the temperature is steeper 
(Figure [2j), as can be seen in panel (a) from Figure [U When the reflected wave reaches the 
cs = va layer it undergoes a secondary transformation, generating a new fast acoustic mode 
and a slow magnetic mode visible below the solid line in Figures [T]i and [U, respectively. 

The fast magnetic mode in the low-/? region is reflected back down due to the gradients 
of the Alfven speed (Figure Hk). If we neglect the contribution of the sound speed to the fast 
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wave speed, the reflection height is given by the layer where the wave frequency uj and the 
horizontal wave number k x are related by uj = VAk X) and it is represented in the figure by a 
dashed line. Around the center of the sunspot the reflection is completed and at the height 
where panel (b) is obtained there is no fast mode in that region. Farther from the axis of 
the sunspot the transformation layer is located at a higher height, and the insufficient height 
of the top boundary of the computational domain avoids the complete reflection of the fast 
mode. 

The Alfven wave can be seen in panels (c) and (/). Along y — Mm, which corresponds 
to (f) = 0°, there is no wave power in the Alfven mode. At this position the magnetic field is 
contained in the xz plane, since B y = G, making this plane equivalent to a 2D case. Under 
these conditions the Alfven mode is decoupled from the fast and slow magneto-acoustic 
modes, and it is not possible for the incident wave to undergo conversion to the Alfven 
mode. As regions farther from the y = Mm plane are considered the conversion to the 
Alfven mode becomes more efficient. It shows a highest amplitude at around 6 = 45°. 
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Fig. 3. — : Logic of the amplitude ratio R = Bi/^/xoPo/vi for the projected directions at 
the upper atmosphere (averages at heights from 0.4 Mm above the c# = va layer up to z = 
1.9 Mm). Left panel: slow acoustic mode (ei ong ); middle panel: fast magnetic mode (e trans )] 
right panel: Alfven mode (e perp ) 



4. Alfven mode polarization relations 



Before evaluating the contribution of the different wave modes to the energy in the 
upper part of the atmosphere, we have checked the validity of the projections described in 
the previous section, in order to ensur e that they provide an ac curate decoupling of the 
slow, fast, and Alfven waves. Following iKhomenko fc Callyl ( 120121 ). we have calculated the 
ratio R = Bi/-^//xoPo/ v i- Since the kinetic and magnetic energy for an Alfven wave is in 
equipartition, one would expect the ratio R to be equal to one for that waves. 



Figure [3] shows the ratio R for the three projected components of the velocity and 
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magnetic field. In each case, Bi and vi pairs were obtained from the decomposition in the 
directions defined by Equations [1] and El The ratio is evaluated at heights from 400 km above 
the transformation layer to the top boundary, and it is averaged in time for the stationary 
stage of the simulations. It reveals the different nature of the three projections. For the 
ei on g direction at 9 < 50°, the ratio R present small values around 10 -2 , indicating that this 
component is dominated by the velocity variations rather than magnetic, and confirming 
that this projection contains the slow acoustic waves. However, at higher inclinations the 
magnetic variations associated with velocity variations are bigger. An opposite behavior is 
found for the e trans component. In this case, in the regions where 9 < 50° the ratio R is 
big, around 10 2 — 10 3 , and it becomes smaller at higher inclinations. These plots indicate 
that the projections ei ong and e trans provide a good estimation of the slow and fast waves, 
respectively, in the upper atmosphere when the inclination of the magnetic field is below 
50°. At higher inclinations both wave modes are mixed up, since the projection directions 
are asymptotic, valid strictly only where the Alfven speed is much higher than the sound 
speed. At these locations this criterion is not fulfilled. 

The right panel of Figure [3] illustrates the ratio R for the e perp component. It shows 
a ratio around 10° for all the atmosphere, indicating that for this projections the velocity 
and magnetic perturbations are in equipartition, which confirm the Alfvenic nature of the 
oscillations in this characteristic direction. 
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Fig. 4. — : Phase shift between the variations of v perp and B perp in the upper atmosphere 
(averages at heights from 0.4 Mm above the c# = va layer up to z = 1.9 Mm) as a function 
of inclination and azimuth (top left panel) and X and Y (top right panel). Bottom panel 
shows the same quantity, but as a function of height and horizontal distance at y — 7.5 Mm. 
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Figure H] shows the phase difference between the perturbations in the velocity and mag- 
netic field in the direction e perp) which correspond to the Alfven wave. For a pure Alfven 
mode a phase shift of 180° indicates upward propagation, while negative phase shifts corre- 
spond to downward propagation. Most of the right part of the sunspot, that is, for azimuths 
between 0° and 90°, presents a phase shift around 180°. At these locations the Alfven waves 
come from the conversion of the upward propagating fast acoustic mode introduced in the 
high-/3 region, and they keep their upward propagation to higher layers. However, at some 
positions in the left side of the sunspot (with (j) between 90° and 180°) the negative sign 
of the phase shift indicates downward propagating waves. The direction of the propagation 
of the Alfven waves in this s imulation shows a perfect agreement with the one obtained 
by iKhomenko fc Callyl (120121 ). Since the efficiency of the conversion to the Alfven mode 
is enhanced with the alignment of the direction of propagation and magnetic field, in the 
right part of the sunspot the upward propagating fast waves couple to upward Alfven waves. 
On the other hand, in the left part of the sunspot the most efficie nt conversion to Alfve n 
waves occurs for t he refracted downward pr opagating fast waves ( jCally fc Hansen! 1201 lh . 
See Figure 1 from IKhomenko fc Callyl (120121 ) for a schematic representation of these mode 
transformations . 



5. Energy fluxes 

The acoustic and magnetic energy fluxes were calculated using Equations [H and El 
respectively. Figures [5] and [6] show the time- average results, including all the time steps after 
the stationary regime is achieved. The former one corresponds to a xz cut at y — 6 Mm. 
The magnetic flux is only plotted above the c$ = va layer, where the distinction of the three 
different modes using the projections described in Equations [TH3] is meaningful. The later 
figure shows a xy plot at z — 1.65 Mm. 

Most of the domain presents a positive acoustic flux, which corresponds to the upward 
propagating slow mode in the low-/3 region. The highest contribution of this mode is located 
at the right side of the sunspot, especially at moderate inclinations where the amplitude of the 
viong projections was higher in Figure [TJ However, at the left part of the sunspot a negative 
flux is obtained. This downward propagating slow mode flux is particularly large near the 
axis of the sunspot and represents the slow waves reflected by the temperature gradient, as 
discussed in the previous section. The 5 mHz frequency fast acoustic wave which propagates 
upward in the high-/? region reaches the conversion layer c# = Va before the high value of 
the cut-off frequency avoids its propagation. Note that near the center of the sunspot the 
cut-off frequency presents its highest value of 5.7 mHz at the temperature minimum, which 
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Fig. 5. — : Vertical component of the energy fluxes at y — 6 Mm. Top panel: acoustic flux; 
middle panel: magnetic flux; bottom panel: magnetic flux due to Alfven waves. The units 
of the color coding are 10 6 erg cm _2 s _1 . Positive fluxes mean energy propagating upward 
and negative fluxes downward. Horizontal solid line is the height where c# = va\ horizontal 
dashed line is the fast mode reflection level. Magnetic field lines are inclined black lines. 



is located 375 km above the c# = va layer. Just above this layer, the recently converted 
slow acoustic mode becomes an evanescent wave. This situation differs from what happen 
at locations far from the axis of the sunspot, where the field inclination reduces the effective 
cut-off and allows the slow mode waves to scape to the upper atmosphere. At a height located 
a few hundreds kilometers higher than the temperature minimum, the cut-off frequency of 
the atmosphere is again below 5 mHz due to the chromospheric increase of the temperature, 
and the slow mode near the axis can propagate again, allowing the downward propagation 
of the waves reflected by the rise of the temperature. 

Above the transformation layer the magnetic energy flux includes the fast and Alfven 
modes. The middle panel of Figure [5] shows that the magnetic flux is positive at most heights, 
while in the middle panel of Figure [6] can be seen how it increases in the Y direction toward 
the periphery of the sunspot. On the other hand, near the X boundaries of the model a 
negative magnetic flux appears, larger than the highest positive magnetic flux. This negative 
flux must be produced by the reflected fast wave. 
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Fig. 6. — : Vertical component of the energy fluxes at z — 1.65 Mm. Left panel: acoustic 
flux; middle panel: magnetic flux; right panel: magnetic flux due to Alfven waves. The units 
of the color coding are 10 6 erg cm _2 s _1 . Positive fluxes mean energy propagating upward 
and negative fluxes downward. 

For a more detailed analysis of the contribution of the different wave modes, the magnetic 
flux of the Alfven wave was deattached from the total magnetic flux by recalculating the 
magnetic flux but using the projections of the velocity and magnetic field in the e trans 
direction in Equation [51 The result is plotted in the bottom and right panels of Figures 
[5] and El respectively. The vertical cut reveals that at y — 6 Mm the vertical Alfven flux 
is positive in the righ part and negative in the left part, that is, at the right part of the 
domain the Alfven wave propagate upward and at the left part they propagate downward. 
This result agrees with the one obtained in the previous section from the phase shift between 
v perp and B perp . The xy cut shows that the Alfven energy flux increases at locations with 
higher inclinations, far from the axis of the sunspot, only for those regions of the atmosphere 
whose azimuth (f) is different from 0°. Thus, the highest Alfven energy flux is obtained near 
the Y boundaries. Note that it is comparable to the highest acoustic energy flux. The 
negative flux at the left side of the sunspot shown in the bottom panel of Figure [5] is hardly 
visible in the right panel of Figure [6] because of the different scale used in the plots. 

A comparison between middle and right panels of Figure El reveals that the negative 
magnetic flux from the X boundaries is missing in the Alfven energy flux, so it corresponds 
to the reflected fast mode, as previously stated. At the positions where the Alfven energy 
flux is larger the total magnetic flux shows a lower positive value, meaning that at these 
locations the fast wave must also contribute with a negative flux. 

6. Energy of the three wave modes 



The wave energy of the slow, fast, and Alfven modes was calculated from Equations EE 
using the corresponding projected velocities. To be consistent with the work by lKhomenko fc Cally 
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( 120121 ) and provide a direct comparison with their 2.5D simulations, time averaged energies 
were obtained at heights from 400 km above the c# = va layer up to the upper boundary of 
the simulation box in the stationary stage of the simulations. 

Figure [7] illustrates the results for all the horizontal positions from the computational 
domain. The left panel shows that there is a prominent maximum of the slow wave energy at 
the right part of the sunspot, not far from the axis. The highest slow mode energy apears at 
around x — 10 Mm and y = Mm. The main interest of this measurement is to quantify the 
amount of energy which is converted at the c# = Va from the incident fast acoustic mode to 
the outgoing upward propagating slow acoustic mode. For this reason, in the left panel the 
regions with negative acoustic flux in Figure [6] have been masked, since the slow wave energy 
in that locations does not comes directly from the conversion layer but from the reflection 
due to the temperature gradient. 

The wave energies have been plotted as function of inclination and azimuth of the 
sunspot magnetic field lines at the c orresponding horiz o ntal l ocations in Figure El This for- 
mat a llows a direct c omparison with |Callv fc Goossensl ( 120081 ) Figure 2, iKhomenko fc Cally 
( 2011 ) Figure 4, and Khomenko fc Cally ( 20121 ) Figure 5. In the left panel the angles corre- 
sponding to the positions with negative acoustic flux has also been masked. It shows that 
the maximum slow wave energy is obtained when the direction of the wave incidence forms 
an angle (f) = 0° with the magnetic fiel d, for an inclination arou nd 6 = 25°. This result 
is consistent with the 2D analysis from ISchunker fc Callyl ( 120061 ) , where they shown that 
around this field inclination the attack angle, i.e., the angle between the wavevector and the 
magnetic field at the conversion layer, is small and produce an enhacement of the conversion 
from the fast acoustic to the slow acoustic wave. The region where the conversion to the 
slow mode is efficient extends to higher azimuths, although it decreases significantly for (j) 
higher than 60°. The distribution of the slo w mode energy as a func tion of 9 and (j) shows 
an excellent agreement with Figure 2 from ICally fc Goossensl ( 120081 ) . The only difference 
appears at higher azimuths, which correspond to the left part of the sunspot s, where the 
reflection of the slow mode produced by the temperature increase complicates the evaluation 
of the upward propagating slow mode energy in this simulation of conversion in a realistic 
sunspot model. 

The Alfven wave energy increases toward the periphery of the sunspot, except along 
y = Mm (right panel of Figure [7j). Its maximum value is obtained very close to the 
boundary of the computational domain around x = Mm. Its representation as a function 
of inclination and azimuth (right panel of Figure [8j) reveals that the conversion to the Alfven 
mode is efficient at inclinations above 40° and azimuths between 50° and 120°. At (j) ~ 50°, 
which corresponds to the right side of the sunspot, the Alfven energy extends to higher 
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inclinations than the left part of the sunspot. This is also visible in Figured The distribution 
of the Alfven e nergy with 9 and <j> i s rem arkably similar to the one found previously from the 
3D analy sis bvlCallv fe Goossensl ( 120081 ) in homogeneous fields. However, a few differences 
arise. In ICally fe Goossensl ( 120081 ) the highest conversion to Alfven waves appears at 9 « 
40°, while in our simulations the maximum energy is shifted toward higher inclinations. A 
similar shift in the position of the maximum was obtained from the 2.5D simulations of 



Khomenko fc Callyl ( 120111 ) in homogeneous magnetic fields. Since the Alfvenic energy peaks 
so close to the boundary of our simulations, it is hard to say if the maximum conversion occurs 
at 9 ~ 47°, as shown by the plot, or at even higher inclinations. Moreover, the maximum 
energy appears in the left part of the sunspot, corresponding to (j) ~ 100°. This differs from 
the analysis of homogeneous magnetic fields , whe re the maximum is loc ated at (j) = 60° 
dCallv fc Goossensl I2OO8I : khomenko fc Callvl boilh or at 6 = 80° - 90° dCallv fc Hansen 



2011 ). In the 2.5D simulations in a sunspot model from Khomenko fc Cally ( 2012h the 
energy of the Alfven wave peaks at azimuths between 70° and 100°. 

Middle panel of Figure [8] shows some fast wave energy at high inclinations. 
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Fig. 7. — : Wave energies at the top of the atmosphere (averages at heights from 0.4 Mm 
above the c# = va layer up to z = 1.9 Mm) for the three projected velocity components as a 
function of the horizontal locations. Left panel: acoustic flux; middle panel: magnetic flux; 
right panel: magnetic flux due to Alfven waves. The units of the color coding are 10 6 erg 
cm -2 s -1 . In the left panel the regions with negative acoustic vertical flux have been masked. 



Discussion and conclusions 



In this paper we have studied the conversion from fast to Alfven modes in a realistic 
sunspot-like atmosphere by means of 3D numerical simulat i ons. This study provid es a direct 
comparison with analy tical works (Callv fc Goossensl 120081 : ICally fc Hansenll201ll ) as well as 



numerical simulations ( IKhomenko fc Cally 



2011 



201211 



With regards to the fast-to-slow conversion, the efficiency of the transformation peaks at 
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Fig. 8. — : Wave energies at the top of the atmosphere (averages at heights from 0.4 Mm 
above the c# = va layer up to z = 1.9 Mm) for the three projected velocity components as a 
function of inclination and azimuth of the sunspot magnetic field lines at the corresponding 
horizontal locations. Left panel: acoustic flux; middle panel: magnetic flux; right panel: 
magnetic flux due to Alfven waves. The units of the color coding are 10 6 erg cm -2 s -1 . In 
the left panel the regions with negative acoustic vertical flux have been masked. 



inclinations around 25° and low azimuths, and i t decreases at higher azimuths. This result 
show s a good agreement with previous works (jSchunker fc Callyl 120061 : ICally fc Goossens 
20081 ) in homogeneous magnetic field configurations. However, the sunspot-like structure of 
this simulation produce some differences in the wave modes of the upper atmosphere. In the 
central and left part of the sunspot the simulation shows downward propagating slow waves, 
which are reflected due to the chromos pheric increase of the tem perature. A similar result 
was found in the 2.5D simulations from iKhomenko fc Callyl ( 120121 ). although in our case the 
downward flux near the axis i s enhanced because o f the higher increase of the temperature at 
the center of the sunspot. In lKhomenko fc Callyl (120121 ) simulations the closest plane to the 
axis of the sunspot is located at a distance of 7.5 Mm, where the increase of the temperature 
in the transition region is lower (Figure [2]). 

The simulation presented in this paper supposes a step forward in the study of the 
conver sion to the Alfven mode, since it accounts for full 3D. The work by IKhomenko fc Cally 
( 120 12l ) was limited to 2.5D, and only a few vertical slices across the sunspot model at several 
distances from the axis were analyzed. Since the fast-to- Alfven conversion depends strongly 
on the angle between the direction of wave propagation and the magnetic field, a full 3D 
simulation is required to retrieve the complete description. The simulation shows that the 
highest Alfvenic wave energy is obtained at high inclinations and for azimuths between 50° 
and 120°. The efficiency of the conversion obtaine d for the sunspot model shows a similar 
pattern to the homogeneous magnetic fields from ICally fc Goossensl ( 120081 ) , although the 
maximum of the magnetic energy of Alfven waves is shifted toward more inclined fields. 



A comparison of the total energy in the upper atmosphere of the upward propagating 
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slow acoustic energy and the corresponding to the Alfven mode reveals that the former one 
is around two times higher. This result remarks the important role of the Alfven energy 
at the chromosphere. At 9 « 25° the acoustic energy is much higher than the Alfven 
energy, especially at low azimuths where its conversion is more efficient. However, at higher 
inclinations the Alfven energy becomes dominant and at 9 = 47° and (f) around 90° it is 
more than 7 times higher than the acoustic energy. The horizontal size of the computational 
domain prevent us to measure the energy at higher inclinations for these azimuths, but from 
the tendency shown in the computed region one would expect a strongest dominance of the 
Alfven energy at higher inclinations. Moreover, the vertical limitation of the domain also 
reduces the energy in the Alfven mode. At the 5 mHz fre quency of this simulat ion the fast- 



to- Alfven conversion region spans over 20 scale heights (jCally fc Hansen! 120 111 ). Since our 



computational box has a much lower size, only a small fraction of the conversion is complete d 



in our domain. The same limitation affects the simulations from iKhomenko fc Callyl ( 120121 ). 
However, in our case we have found a higher energy in the Alfven wave due to the higher 
inclinations covered by our computational domain. 

Despite of the dependence of the conversion with the azimuth, due to the stochastic 
direction of the wave propagation in the Sun, one would expect to find an efficient conversion 
to the Alfven wave at all the locations surrounding a sunspot with a certain inclination. Near 
the umbra-penumbra boundary, the waves whose attack angle is below 90° will generate 
upward propagating Alfven waves, while the incident waves which propagate in the opposite 
direction will produce downward Alfven waves. In this way, in these regions upgoing and 
downgoing waves will coexist. On the other hand, at the higher inclinations of the penumbra, 
the energy flux of the Alfven mode is only composed by upward propagating waves for all 
the directions of incidence of the fast acoustic wave where the conversion to the Alfven wave 
is efficient, as shown by the right panel from Figure [6] at 9 « 45°. It differs from the situation 
at 9 « 30°, where at some regions of the sunspot there is negative flux. In the case of solar 
observations, where the incident waves propagate in all directions, the Alfven waves in an 
annular region surrounding the sunspot with higher inclinations (around 9 ~ 45°) will only 
consist on upward propagating waves. 

The small time step imposed by the high chromospheric Alfven speed has been an 
important concern in the development of the numerical simulations presented in this paper. 
Some considerations have been taken into account in the numerical configuration regarding 
this issue. Firstly, a magnetic field strength of 900 G was adopted at the axis of the sunspot 
at the photosphere. This value is well below what might be expected in a mature sunspot. 
A more realistic field strength would lower the height of the c# = va and the fast mode 
reflection layers. This would allow the fast-to- Alfven conversion to be produced in a bigger 
region of the computational domain, and might generate an even higher Alfven energy flux 
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at the top of the computational d omain. Secondly, this t op boundary has been set at the 



chromosphere. The recent work by lHansen fc Callyl (120121 ) has shown the interesting effects 



of the chromosphere-corona transition region on fast-to- Alfven conversion. They found that 



the reflection of the Alfven waves at the transition region (lUchida fc Sakurailll975l ) is sensitive 
to the distance between the fast reflection point and the transition region. The Alfven flux 
that can reach the corona is increased when this distance is small. In the present work the 
computational box does not include the transition region, since it would produce an even 
smaller time step, compromising the realization of this simulation using a reasonable amount 
of computational time. It would be interesting to extend the simulations to the corona and 
adopt a realistic magnetic field strength, but it is out of the scope of this paper and shall be 
accounted for in forthcoming works. 
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for their examination and suggestions to an early version of this paper. This work used 
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